Branching, Capping, and Severing in Dynamic Actin Structures 
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Branched actin networks at the leading edge of a crawling cell evolve via protein-regulated pro- 
cesses such as polymerization, depolymerization, capping, branching, and severing. A formulation 
of these processes is presented and analyzed to study steady-state network morphology. In bulk, 
we identify several scaling regimes in severing and branching protein concentrations and find that 
the coupling between severing and branching is optimally exploited for conditions in vivo. Near the 
leading edge, we find qualitative agreement with the in vivo morphology. 
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When a cell crawls, it must reorganize the cytoplas- 
mic network of biopolymers that controls its shape. 
The shape and motility of the leading edge (the lamel- 
lipodium) of a crawling cell are determined primarily by 
a dynamic network of actin filaments (F-actin) [l[ . These 
filaments are living polymers made up of monomeric G- 
actin, and have a definite polarity such that monomers 
tend to be added to the plus (or "barbed") end, and 
tend to fall off at the minus (or "pointed" ) end. A host 
of regulatory proteins concentrate plus ends at the lead- 
ing edge by controlling the morphology of the actin fil- 
ament network 2]: capping proteins bind to the plus 
end of filaments, preventing further growth; severing pro- 
teins, such as cofilin, bind to filaments and break them in 
two, enhancing depolymerization; and the protein com- 
plex Arp2/3 nucleates branches from existing filaments, 
creating new plus ends at each branch tip and protect- 
ing the bound minus ends of branches from depolymer- 
ization. These proteins and their associated processes 
constitute the dendritic nuclcation model 2]. This bi- 
ological model is supported by experiments that show 
that severing, capping and branching are all crucial to 
cell motility [a, Hi, HQ- 

In order to gain insight into how the cell controls 
its structural dynamics it is imperative to understand 
quantitatively how actin morphology can be controlled 
by the concentrations of these regulatory proteins. To 
date, theoretical studies of these kinetic processes have 
either omitted one or more of these critical processes 
@, S EH EH, Eli or have been restricted to describing 



just the overall amount of polymerization as a function 
of time Other treatments of actin-based motility 

have focused on the interplay of force generation with 
motility 14 , EB| ; and treat the branched actin network 



either as an elastic continuum or as a collection of un- 
coupled filaments. In this paper, we adopt a complemen- 
tary approach that not only captures all of the critical 
features of the dendritic nucleation model but also allows 
the first theoretical investigation of the morphology of the 
branched network. We confirm biological and biochemi- 
cal understanding of the proteins involved by obtaining 
quantitative agreement with bulk in vitro experiments. 



We extend the approach to explore the interplay between 
actin morphology and motility near the leading edge of a 
crawling cell. The agreement of our theory with electron 
microscopy images of the leading edge of a crawling cell 
provides strong evidence that the dendritic nucleation 
model captures the key players in this important form of 
cell motility. 

We introduce a mean-field rate equation formulation 
for actin structures that undergo polymerization, de- 
polymerization, capping, branching, and severing. We 
capture key morphological information by retaining the 
entire length distributions of filaments with free minus 
ends and of branches. Our formulation is constructed for 
polar actin assembly with branching (Arp2/3), severing 
(cofilin) and capping agents, whose action is described by 
rate constants. Initially, capping is taken into account by 
assigning a probability for a filament to be capped; for 
reversible capping, this renormalizes the polymerization 
rate. Thus, k + is the effective polymerization rate con- 
stant. The growth rate at the barbed end of a filament is 
k + p m , where p rn is the free monomer concentration. The 
depolymerization rate is denoted by . Nuclcation and 
dissociation of filament seeds are described by rate con- 
stants k n and kd [IB| respectively. Nucleation of branches 
is modeled by the rate constant k arp = k^ rp pA rp 2/3Pm 
[liL E3], where /OAr P 2/3 is the concentration of Arp2/3. 
Eventually, Arp2/3 dissociates and branches detach at 
a rate kdr [13] • Severing occurs at a rate k s = k° s p s , 
where p s is the severing protein concentration [lq . The 
nucleotide state of the polymerized actin controls the 
rates of both branching (preferred at ATP-actin sites[l9|) 
and severing (preferred at ADP-actin sites fiojl). We 
therefore consider hydrolysis and phosphate release oc- 
curring in sequence to convert newly polymerized ATP- 
actin to ADP-actin at a rate k pr [2l|. The probability 
that a monomer at distance L from the barbed end of 
a filament is an ADP-actin monomer is then given by 
p(L) = 1 — exp(— L/l c ) [9], where l c — k + p m /k pr . We 
allow branching only on ATP-actin monomers and sev- 
ering only on ADP-actin monomers in filaments. 
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Unless specified, k + = 8.7 pM' ~s 
0.03 s" 1 (estimates are of the order of 
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[23| depending on conditions), k d = 5 x 10~ 4 s _1 minus ends (i.e. minus ends that can depolymerize), 

k n = 3 x 10 _6 /iM _1 s _1 (fitted; see caption of Fig. [T]), Pu(L), (2) the density of branches of length L (which have 
k® — 4.7 x 10~ 4 pM~ 3 s~ 1 (fitted and consistent with bound minus ends that cannot depolymerize), p b (L), and 



M) k dr = 0.0018 s"LJp, UM, = 4 X 10 " 5 mM"^" 1 (3) the density of monomers, p m . In the mean-field bulk 

(estimated from [3, ij) and k pr = 0.002 s" 1 (2^. To case where all rate constants and densities have no spatial 

capture the morphology we consider (1) the density of dependence, we have: 
actin filaments of length L (in monomer units) with free 



L-l 



Pu(2) 
Pfc(2) 



p'u(L) = -k + p m (p u (L) - p u (L - 1)) + k-(p u (L + 1) - Pu{L)) + k dr p b (L) - k s ^ p(L')p u (L) 

L' = l 

oo oo 

+ J2 k sP (L) Pb (L') + J2 k s (p(L)+p(L'-L))p u (L') 

L'=L+2 L'=L+l 

L-2 oo 

Pb(L) = ~k + p m (p b (L) - p b (L - 1)) - k dr p b (L) - k s ^2 P( L ')Pb(L) + k s p(L' - L)p b (L') 

L' = l L'=L+l 
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L=2 L' = l 



L'=3 



withp(i') = p{2) +p(L' — 2). A fifth equation is the con- 
servation of total number of monomers. These equations 

are similar to those in Refs. [nimHE!, 

though we 

distinguish between unbranched and branched filaments, 
which is necessary to quantitatively describe morphology. 

We first study the steady state in absence of severing 
and nucleotide dependence. The length distributions are 



p u {L)=Ae- L / K +Bc 
Pb (L) = Ce~ L ^ 
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where A „ = 1/ log 
1/ log kdr+k + Pm ~ k + p " 



A,B 



fc_- 

and 



(5) 
(6) 

A h = 



k + p„ 

C are length- 



fc+Pm k dr 

independent and depend on the rate constants and total 
actin concentration. The approximations are valid for 
L » 1, which is generally the case in vivo, and imply 
that k- « k + p m and Ah = k + p m /k dr w k-/k dr . 

We performed Brownian dynamics simulations to de- 
termine the accuracy of Eqns.(5,6) (26[. We have 
used rates somewhat different from physiological ones to 
achieve equilibration within reasonable time. Fig [IJa) 
shows that the simulations and theory yield length dis- 
tributions in good agreement with each other with no 
adjustable parameters at these concentrations. The 
predicted maximum in p u (L), which occurs because 
branches can fall off, is also observed in the simulations. 

We have also compared our calculations to in vitro ex- 
periments by Blanchoin, et al. 27j. We neglect sev- 



ering, depolymerization and nucleotide-state dependence 
because those experiments did not contain severing pro- 



teins and filaments were stabilized by phalloidin. Our 
results for the barbed end concentration are plotted as 
a solid curve as a function of Arp2/3 concentration in 
Fig. [Hb), in good agreement with the experiment. 

When severing and monomer nucleotide state are in- 
cluded, analytical progress can be made by using "global" 
conservation principles at steady state to derive simple, 
approximate expressions for morphological properties. 
In steady state, the average total number of branches 
must be conserved. Equating the rate of production 
and destruction of branches gives k dr J^^°_ 2 Pb(L) = 

k arpPArp2/3Pli Ei=2 (Ei'=l ( 1 -p( L ') ) ) (Pu {L)+Pb (L) ) M 



karpPArp2/3f(-jT^ 



\2 tot 
I Prt 



where p£ is the total actin con- 
centration and / denotes the fraction of polymerized F- 
actin monomers that are ATP-bound and hence capable 
of supporting branching and immune to severing. 

Similarly, the average number of free filaments is 
fixed such that k s (l - /) J2T=2 L(pu(L) + Pb{L)) + 
k dr ^2iT, = 2 Pb(L) = k d p u (2). If we assume that 
p{L) can be described by one characteristic length, 
L u , implying p u (2) 

l k d 



'ks(l-f)+k arpPArp2/3 f(^ 



p^/L u , then L u = 
In steady state, conserva- 



tion of ATP F-actin implies / « (^pr-jr 14 + l) such that 



Lu - ( 
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where k' arp = kP arp p Arp2 /z- Comparing the two terms 
in the denominator, we see that branching does not sig- 
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FIG. 1: (a) Length distribution of filaments and of branches 
(in monomer units), as predicted by our mean field the- 
ory (lines) and as calculated from Brownian dynamics sim- 
ulations (points), with k+ = 15 /xM~ 1 s~ 1 , k- = fc<j r 
I UN s ,k a rp/k+ = 0.015. (b) The barbed end concentration 
as a function of Pai- P 2/3) as predicted by mean field theory 
(lines) and measured by experiment (points) [27[. Here, k n 
and k® rp were estimated by fitting to average filament length 
and branching percentage data from [2tJ (not shown). Inset: 
The predicted steady-state value of the fraction of ATP-bound 
actin monomers in filaments, as a function of L u for various 



nificantly affect the average length L u until the Arp2/3 
concentration, pArp2/3, is at least four orders of magni- 
tude larger than the cofilin concentration, p s . In vitro 
motility experiments 0] typically operate in the regime 
where /?Arp2/3 < Ps- Then Eq. [7] shows that for large 
k s , L u is short and scales as L u ~ (kdk-/(k s k pr ))i . 
For lower k s , (but k s » k' arp (^) 2 ), L u is larger and 

shows a more sensitive dependence on k s : L u ~ (kd/k s )i. 
The crossover between these two regimes occurs when 
L u ~ 10. On the other hand, when severing is negligible 
and k pr is small such that side-branching dominates, then 
L u oc 1/ ■ s /k' arp . For small k s and large phosphate release 
rates k pr , where end-branching dominates, L u oc l/k' arp . 
These scaling predictions are experimentally testable. 

The density of branched filaments is also constant in 
steady state. As k pr increases, branching occurs only 
near the barbed end of filaments. Thus, the branch den- 
sity will depend not on the number density of monomers 



in filaments, but on the number density of filaments, 
which is sensitive to severing and is given by Pb oc 
\fk~ s k' arp . On the other hand, if there is no dependence 
on nucleotide state, the branch density depends on the 
total F-actin concentration, which depends on k' arp , but 
not sensitively on k s . 

These results suggest that the dependence of branch 
density on the nucleotide state of filament monomers al- 
lows the cell to control morphology by varying the con- 
centration of severing protein [3] • For in vivo values of 
k pr , the system switches from a side-branching regime 
(/ = 1) to an end-branching regime (/ = 0) as the av- 
erage lengths increase from tens to several hundred sub- 
units, the range typically observed in vivo. In Fig.QTb), 
we plot the predicted value of / as a function of average 
filament length for different values of the phosphate re- 
lease rate. For the in vivo value of the phosphate release 
rate and the in vivo range of filament lengths (ranging 
from tens to several hundreds of monomers), the range 
of / that is covered is almost maximal. This suggests 
that in vivo, the branching morphology (controlled by /) 
is maximally sensitive to severing proteins. Branching 
and severing have been viewed as antagonistic, because 
branching promotes filament growth and severing pro- 
motes depolymerization. However, our results show that 
all processes combine to give maximal cooperativity be- 
tween branching and severing. 

We now consider the coupling of morphology to motil- 
ity by confining the system between two hard walls both 
moving with velocity v. The "front" wall models the 
leading edge of the cell and while the "back" wall de- 
notes the the back edge of the lamellipodium. Because 
capping is necessary to motility Q, we extend Eqns. (1- 
4) to explicitly include a capping rate k c and uncapping 
rate k uc by defining four populations of filaments: capped 
and uncapped branched and unbranched filaments. Now 
p u {L, t, z) denotes the density of capped or uncapped fil- 
aments of length L whose barbed ends are at distance 
z from the front wall. Experiments find that filaments 
are typically at an angle of 9 = 35° with respect to the 
leading edge (28|, so we assume that all filaments are at 
that angle with the wall. Because Arp2/3 is activated 
at the leading edge and diffuses away from it, we model 
the Arp2 /3 concentration profile as an exponential decay 
from the front with a decay length z arp [29]. Similarly, 
we use a severing concentration profile that is an inverted 
exponential with a rise length of z s to model the effect 
of nucleotide dependence on the severing efficiency. Fi- 
nally, the amount of free monomers is determined by the 
conservation of actin, and actin that flows beyond the 
back wall is recycled in as free monomers; our results 
are for systems sufficiently large such that actin reaching 
the back wall is small compared to actin within the sys- 
tem. We use an exponential profile with a decay length 
of D/vt, consistent with the diffusion equation, where D 
is the diffusion constant of free monomers and Vt is the 
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FIG. 2: (a) L vs. distance from leading edge, z, for differ- 
ent cofilin concentrations, p s . Here, v — 0.1/im/min, t>t ~ 
25 /im/s, D — 5 /im 2 /s, z arp = 0.1 /im, z s = 2 /^m, k c = 1 s _1 , 
fc„ c = 0.01 s" 1 with 1 /iM Arp2/3 and 100 fiM actin (at the 
upper end of the physiological range), (b) Distance between 
branches vs. z for the same p a as in (a). Inset: Lamellipodium 
length vs. p s . 

transport (bias) velocity of free monomers towards the 
surface. We choose v t = 250 v, where v is the speed of 
the moving surface, consistent with results of Zicha et 
al. [3(3]. Not all values of v t yield reasonable morpholo- 
gies. If, for example, v t = 500 v, the system tends toward 
a state with no filaments. 

In steady-state, we find that branching, severing, and 
capping are all needed to obtain a morphology that is 
consistent with experiments. Branching in the front is 
needed to help the system "keep up" with the wall by nu- 
cleating new filaments. Capping towards the front chan- 
nels new filament growth into shorter branches and thus 
increases filament density. Finally, severing is needed to 
replenish the free monomer supply. Fig. [2] shows the 
steady state morphology for the system with a moving 
surface. In Fig. EJa), we plot the average length of fila- 
ments (branched and unbranched) as a function of dis- 
tance z from the front wall. The average length, L, is 
short near the surface, and increases to approximately 
100 nm at a distance of roughly 0.5 //m from the surface. 
As the severing protein concentration p s increases, L de- 
creases with increasing z, eventually reaching zero when 
all filaments are depolymerized away. The distance at 
which L reaches zero provides an upper bound on the 



lamellipodium length, i .e . the distance the lamel- 
lipodium extends into the cell from the leading edge. 
This length is plotted in the inset to Fig. [2|b) as a func- 
tion of p s . For tens of micromolar concentrations, the 
lamellipodium length is tens of microns, consistent with 
experimental observations 3l|. The predicted length in- 
creases with actin concentration. 

Fig. [2jb) shows that the average distance between 
branches monotonically increases with z. The higher p s , 
the greater the increase in the branch distance with z. 
In short, Fig. [2] implies that the morphology consists of 
short branched filaments within the first micron of the 
leading edge, with longer, less branchy filaments further 
away. These observations are consistent with electron 
microscopyimages of the lamellipodium in crawling ker- 
atocytes [31|, suggesting that the dendritic nucleation 
model captures the minimal set of proteins involved. A 
thorough, quantitative experimental test of our results 
near a moving surface should provide a stringent check 
on the validity of the dendritic nucleation model for actin 
assembly near the leading edge of a crawling cell. 
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NSF-DMR-0520020. 
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